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ABSTRACT 

We carry out a detailed study of the orbital dynamics and structural evolution of over 6000 subhalos in the 
Via Lactea II simulation, from infall to present. By analyzing subhalos with masses down to m = 4 x 10^ M0, 
we find that lower mass subhalos, which are not strongly affected by dynamical friction, exhibit behaviors qual¬ 
itatively different from those found previously for more massive ones. Furthermore, there is a clear trend of 
subhalos that fell into the host earlier being less concentrated. We show that the concentration at infall charac¬ 
terizes various aspects of subhalo evolution. In particular, tidal effects truncate the growth of less concentrated 
subhalos at larger distances from the host; subhalos with smaller concentrations have larger infall radii. The 
concentration at infall is further shown to be a determining factor for the subsequent mass loss of subhalos 
within the host, and also for the evolution of their internal structure in the Vmax “ ^max plane. Our findings raise 
the prospects of using the concentration to predict the tidal evolution of subhalos, which will be useful for 
obtaining analytic models of galaxy formation, as well as for near field cosmology. 


1. INTRODUCTION 


The standard model of cosmological structure formation is 
based upon the notion that the gravitational landscape is domi¬ 
nated by cold dark matter (CDM), which initially collapses on 
small scales and grows hierarchically to larger scales through 
the continued merger and accretion of smaller objects. Many 
of the accreting systems survive to the present epoch as inde¬ 
pendent entities within their host, giving rise to a system of 
nested substructure within the largest objects to have formed 
today. These surviving remnants present a unique opportunity 
to study the seeds of galaxy formation and the chance to probe 
the nature of dark matter on small scales. 

Accomplishing these goals requires the development of an 
accurate and predictive theory for the evolution of substruc¬ 
ture. This remains a difficult task and a fundamental problem 
in the burgeoning area of near field cosmology. Even without 
taking into account the complicating role played by dissipa¬ 
tive effects associated with, e.g., star formation and feedback, 
one must still contend with reconciling the stochasticity in the 
primordial fluctuations that seed the subsructure in any given 
object, on the one hand, and the highly nonlinear gravitational 
dynamics associated with tidal disruption and dynamical fric¬ 
tion, on the other hand. As is often the case, efforts to tackle 
this problem generally fall into one of two cat egories. _ 

First are direct numerical simulations (e.g., [Diemand et al.| 
20081 |Springel et al.||2008| |Garrison-Kimmel et al.||2014|l, 
which attempt to solve the problem ab initio from cosmolog- 


ical initial conditions zoomed on a single host. This class of 
approaches suffers from small number statistics, both in the 
number of individual objects simulated and in the range of 
underlying background cosmologies (warm dark matter, self¬ 
interacting dark matter, broken scale invariance, etc.). In ad¬ 
dition, while simulations have begun to converge on accurate 
solutions for individual systems in a ACDM universe com¬ 
posed only of collisionless dark matter, star formation and 
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feedback must still be treated using heuristic sub-grid ap¬ 
proaches calibrated with empirical data, blurring the line be¬ 
tween theory and observation and complicating the interpre¬ 
tation of simulations. 


The second class involves semi-analytical galaxy formation 

models (e.g. 

Taylor & Babul|2004 

jZentner et al.||2005| |Gan[ 

jet al.|2010[[Jiang & van den Bosch 

2014b| [Pullen et al.|2014|l 


which have a long history in cosmology, and have begun to 
be successfully applied to the local universe. The standard 
approach is to generate a mass accretion history using the ex 

idividuc 


cursion set formalism (e.g.. Bond et al 


Lacey & Cole 

199311 followed by an integration of individual accreting or¬ 


bits from the moment of infall to the present day. Orbital 
parameters at infall are drawn from pro bability distributions 
motivated by numerical simulations (e.g.,|Navarro et al.|1995[ 


;n| 19971 

Wang et 


Ghigna et al.|1998||Benson|2005||Zentner et al. 

2005 Wang et al.|2005[ kliochtar & Burkert|2006[ Jiang et al. 

2008 Wet^l|2011[ Jiang et al.|2015p and the time integration 

contains prescnptions for various nonlinear processes such as 
tidal stripping, tidal heating, and dynamical friction. The util¬ 
ity of this approach is its computational speed, allowing one 
to simulate multiple realizations and cover a broad region of 
model and parameter space compared to what can be achieved 
using expensive cosmological simulations. The drawback is 
that any simplifying assumptions (e.g., symmetries in the host 
potential, omission of substructure interaction) inherent to the 
model may affect the final result in an unknown or unphysical 
manner. 

In this paper, we present a detailed case study of an indi- 
vi dual object simulated at high resolution - the Via Lactea 
II ( [Diemand et al.||?008| VL2) simulation - with the aim of 
making connections relevant to semi-analytic models of sub¬ 
structure evolution. We focus on a self-consistent description 
of the most important physical processes and relationships, 
rather than on direct comparison to specific observations. Our 
goal is to separate the robust quantitative predictions of this 
simulation from those that are unique to the particular back¬ 
ground cosmology and random realization used to generate its 
initial conditions. 

This paper is organized as follows. In Section we de¬ 
scribe our methodology of extracting substructure evolution 
from the public VL2 catalogues. In Sectionj^we present the 
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main results of our work. We begin in §3.1 with a basic de¬ 


scription of the host halo and in ^3.2| statistics of its subhalo 
population, followed in §3.3| with a pr esen tation of orbital 
properties at the time of infall, and in §3.4| with a quantita¬ 
tive assessment of substructure evolution including the physi¬ 
cal processes of tidal mass loss and its dependence on subhalo 
properties, the orbital timescale, changes in the orientation of 
the orbital plane, and the dynamical readjustment of the inter¬ 
nal structure of subhalos as portrayed by their movement in 
the rmax-Vmax plane. In ^we compare the surviving and dis¬ 
rupted subhalo populations of VL2 and investigate how sur¬ 
vivability depends on infall redshift, mass, concentration, and 
orbital parameters. We summarize our conclusions in Section 

El 

2. DATA ANALYSIS 


The VL2 simulation traced the growth of a galactic host 
halo within a high-resolution region sampled with roughly 
one billion particles of mass 4100 Mq. In what follows we 
make use of the main halo catalogue made publicly availably 
by the VL2 team. This catalogue contains evolutionary tracks 
of all 20048 (sub)halos within the simulation box that are re¬ 
solved at z = 0 and for which their peak circular velocity was 
larger than = 4 km s“' at some time during their evolu¬ 
tion. The latter restriction is imposed to discard small halos 
affected by insufficient resolution. 

The catalogue contains a collection of halo properties at 
27 discrete redshifts between 0 < z < 27. These properties 
include: the x, y, and z positions and velocities relative to 
the host halo rest frame; the tidal radius, rdd, and tidal mass, 
iMtid; the maximum of the circular velocity curve, Vmax^ and 
the radius, r^ax, at which this occurs. Empty values occur at 
redshifts when the halo progenitor either did not exist or over¬ 
lapped with a more massive halo. In what follows we consider 
only the redshift range 0 < z < 4.56 for which the host pro¬ 
genitor was consistently identified within the simulation. This 
contains 19 redshift snapshots which we further refine by per¬ 
forming cubic spline interpolations of the above quantities to 
generate a total of 181 discrete sample points equally spaced 
by 68.8 Myr. 

Subhalos are identified in VL2 using the six-dimensional 
p hase-space friends-of-f riends (6DFOF) algorithm described 
in |Diemand et al.| ( |2006| l. Around each (sub)halo the circular 
velocity profile, Vcirc = \/Gm{< r)/r, is computed in spherical 
bins and is fitted with the sum of contributions from an NFW 
profile and a constant density background, pbg. The latter 
component is then subtracted from the (sub)halo density pro¬ 
file and a tidal radius is computed by solving /5sub(i‘tid) = 2pbg, 
corresponding to the tidal radius of an isoth ermal sphere on 
a circu lar orbit within an isothermal host ( |Diemand et al.| 
2007 aj. The tidal mass is assigned mbd = m{< rdd). For suffi¬ 
ciently isolated halos, where the background density is small 
and rtid > r 2 oo (the radius at which the enclosed density is 
200 times the mean matter density), rdd is capped at r 2 oo and 
'Mtid = 'W 200 (Diemand, private communication). 

Hence, the subhalo masses used in this paper are not the re¬ 
sult of an unbinding procedure of dark matter particles. Nev- 
ertheless, it was shown in the Via Factea I (VFl) analysis 
( [Diemand et al.|[2()07b[ ) that this definition of tidal mass in¬ 
deed agrees well with the true bound mass when the subhalo 
is near apocenter, but may significantly underestimate bound 
mass near pericenter. For this reason we generally only report 

^ http://www.ucolick.org/^diemand/vl/data.html 


mass quantities near apocenter and explicitly point out to the 
reader when this is not the case. 

In the following subsections we define concepts and present 
our methodology of investig ating substructure evolution from 
the VF2 data. We begin in |2.1| with the definition of a sub¬ 
halo. In §2.2|we model the internal structure of the host and 


its subhalos via concentration parameters. We define in §2.3 


the redshift, Zinfaih at which a subhalo is said to first infall onto 
the host. In §2.4| we outline our calcul ation s of orbital energy 
and angular momentum and finish in §2.5| with a description 
of our method of tracing subhalo orbits after infall. 

As a matter of convenience, we remove explicit redshift 
dependence in our following notation and remind the reader 
here that all quantities are computed at discrete times. We 
use lower case notation (e.g., m, rm^x, Vmax) when referring 
to subhalos while upper case notation (e.g., M, /?max, Vmax) 
is reserved for the host. The mass of a subhalo is taken to 
be its tidal mass while the mass of the host is taken to be its 
virial mass (see §2.2 1 . We often use /r = m(z)/M{z) to de¬ 
note the instantaneous mass ratio between a subhalo and the 
host. At times we normalize to the present-day host mass in 
which case we define = m(z)/M{Q). In what follows we as¬ 
sume the same cosmology as the VF2 simulation; namely, the 
ACDM parameters (H^, Ha, h, «j, erg) = (0.238, 0.762, 0.73, 
0.951, 0.74) from the WMAP 3-year data release (ISpergel 
jet al.|2007| l. 


2.1. Definition of a subhalo 

We flag an object in the VF2 catalogue as a subhalo if at one 
time during its evolution it passed within the instantaneous 
virial radius of the host. This definition includes subhalos 
that are presently within the virial radius as well as subhalos 
that currently reside outside the virial radius. We refer to the 
latter group as ejected subhalos in the sense that they are now 
removed from the virial boundary of the host. This terminol¬ 
ogy, however, does not impl y that these subhalos are unbound 
from the host, as shown in §3.3.3| 

In late r sections we explore sub halo tidal mass loss. It was 
shown in jKazantzidis et ar] ( |2004| ) that subhalos with too few 
particles within their tidal radius experience artificially large 
tidal mass loss. For this reason we impose a further restriction 
on the VF2 catalogue that only objects with at least 100 parti¬ 
cles in their tidal radius at z = 0 may be considered as subha¬ 
los. This sets a minimum mass resolution of m = 4 x 10^ Mq. 

We find a total of 7569 objects meeting the above criteria. 
5845 (77%) of these currentl y re side within the host virial 
radius of /?vir = 320 kpc (see ^2.2| l while the remaining 1724 
(23%) are currently ejected. For the remain der o f the paper we 
exclude those subhalos whose infall (see ^2.3| ) is determined 
to be Zinfaii > 4.56. This reduces the total population to 6145 
subhalos with 4607 (75%) currently within the virial radius 
and 1538 (25%) ejected. 


2.2. Host and subhalo mass distributions 


It was shown in [Navarro et al.j ( [1997[ ) that dark matter ha¬ 
los obey a universal density profile, named an NFW profile 
after its founders. This has the form p{r) cx where 

X = r/rs and rj is the radius at which d\np/d\nr = -2. The 
virial radius, r^\^, is defined such that the enclosed density is 
A(z) times the critical density, where A(z) is calculated using 
the fitting function to the overden sity of a virialized uniform 
sphere in a flat universe given in [Bryan & Norman (19981. 
An NFW profile is often parameterized by its concentration. 
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Cvir = ?"vir/rj, which describes the degree to which the mass is 
contained within the central region. 

We assume that the density profile of the host follows an 
NEW form. We determine its concentration by finding the 
unique NEW profile for which the mass enclosed within R^ax 
is ^maxi^max/^- involves the implicit solution of 


^(Cvir) = g(Xm) 


A(z) 


H(z)R„ 


n 2 


Vm 


( 1 ) 


where g(x) = f(x)/x^ with /(x) = ln(l +x)-x/(l +x), //(z) is 
the Hubble parameter, and x„, = f?max/G ~ 2.163. Once the 
concentration is obtained, the host halo mass is computed as 

Afvir = ^TTPeril where f^vir “ ^vir^max/ 

The assumption of an NEW profile for the host should be 
valid over the entire redshift range considered here. Eor sub¬ 
halos, however, an NEW profile is only valid up until its mo¬ 
ment o f infall onto the host. It was shown by [Hayashi et al.| 
( |2003| ) that the processes of tidal heating and stripping tend 
to modify the internal stmeture of subhalos away from their 
initial form. Eor this reason, we only use Cvir obtained from 
equation Q for subhalos at their time of infall. Afterwards, 
we define a concentration parameter 


^max 


= 2 


^max 


2 


Ho>'max 


( 2 ) 


which gives the mean density within r^ax in units of the crit¬ 
ical density. Comparing to equation 0 shows that, for any 
given redshift, there exists a monotonic relationship between 

Cvir and Cmax- 


2.3. Definition of infall 

In an idealized description, a subhalo will form distinct 
from its future host, accreting surrounding material and grow¬ 
ing steadily in size. This process will occur until the time at 
which tidal interactions with its host become important. At 
this point, the combined action of dynamical friction and tidal 
stripping will cause the subhalo to lose mass over time. We 
therefore define infall as this turnaround phase in the growth 
history of the subhalo. That is, we define the redshift, Zinfaii, 
of infall onto the host to be the moment in time at which the 
mass of the subhalo is a maximurr|£] 

Recall that masses in VL2 are assigned as the mass con¬ 
tained within the tidal radius. The tidal radius is derived, at 
any moment, by equating the subhalo density profile to twice 
the local background density. Eor sufficiently isolated sub¬ 
halos, where the background density is small, the tidal ra¬ 
dius is capped at r 2 oo- The resultant tidal mass provides a 
good estimation of the true bound mass when the subhalo is 
near apocentre, which is generally the case at infall. As men¬ 
tioned earlier, we consider only the 6145 subhalos for which 
Anfaii < 4.56 since at earlier times the host progenitor is only 
sporadically identified within the VL2 catalogues, preventing 
us from computing orbital properties at those times. 

An alternative convention that is commonly used in the lit¬ 
erature is to define infall as the moment the subhalo passes 

^ This definition does not filter out the possibility that a subhalo may ini¬ 
tially lose mass via tidal interactions wit h a halo other than i ts future host. 
Such group preprocessing was studied in|Wetzel et al.|(|20151 where it was 
found that a significant fraction of subhalos reside within the virial radius of 
a different halo prior to passing through the virial radius of the main host. We 
therefore note the possibility that our zinfaii are biased toward larger values, 
though it is unclear to what magnitude group preprocessing affects premature 
mass loss. 


throug h the virial radius of the host. However, as shown previ¬ 
ously (|Hahn et al.|2009| |Behroozi et al.|2014|l, subhalos gen¬ 
erally undergo strong tidal forces at distances larger than Rvir- 
Eurthermore, the virial radius evolves with redshift through 
its dependence on A(z), meaning that its value will change 
even if the intrinsic mass profile of the host is unchanging. 
The virial radius is therefore not well-suited for defining the 
distance at which a subhalo becomes tidally truncated by the 
host, and can be said to undergo infall in the sense considered 
here. 


2.4. Orbital energy and angular momentum 

We determine the energy and angular momentum by assum¬ 
ing that subhalos evolve as isolated point particles within the 
spherically symmetric NEW profile of the host. In this case, 
the host potential is 

— ^maxVmax ltt( 1 

fiXm) r 

where r = |r| is the radial separation between the subhalo and 
host. In this expression we have taken the zero point of the 
potential to be at infinity. The specific orbital energy of the 
subhalo is 

£■=-vv-l-$(r), (4) 

while its specific orbital angular momentum is 

L = r X V. (5) 


Here v is the physical relative velocity between the subhalo 
and host which includes the sum of peculiar motion and Hub¬ 
ble flow. 

In general, the continued action of dynamical friction will 
steadily drain energy and angular momentum from subhalo 
orbits. The evolution of a subhalo after infall thus depends 
strongly on its initial energy and angular momentum. It is 
therefore important to characterize the infall distributions of 
E and L as inputs in semi-analytic models of subhalo evolu¬ 
tion. A common convention is to normalize these quantities 
in terms of a circular orbit of the same energy. We introduce 
two variables: (i) rj = rcuc/Rm, defined to be the ratio of the 
radius, rcirc, of a circular orbit of the same energy E as the 
subhalo to the virial radius, Ryir, of the host at infall; (ii) the 
circularity, e = L/Lcirc, defined to be the ratio of the subhalo 
angular momentum, L, to the angular momentum, Lciic, of a 
circular orbit of the same energy. To compute rj and e we must 
first evaluate r^jc, which is achieved by numerically solving 
the expression 

ln(l+y) 1 _ 2Ef(x„) 

V 1-i-v X ^ 

y 1-ry •^m^'max 

where y = x„r^i,jR^^. Then Lcirc = where 

Mirdrc) is the mass contained within radius rcirc of the host. 

The definitions of t] and e used here are self-consistent with 
the description of a subhalo orbiting within an isolated NEW 
profile. This does not, however, conform with the standard 
method applied in semi-analytic models of substructure evo¬ 
lution. Instead, it is common to report these quantities at the 
time when the subhalo first crosses through Ryir and to model 
the host potential as a point mass of Myir. In this case, the 
orbital energy is 

1 GMvi, 

E= -V-V--— 


(7) 
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Fig. 1.— Evolution of the viral mass (top panel) and concentration (bot¬ 
tom panel) of the host halo obtained by finding the unique NFW profile that 
matches the values of R^ax and V^ax at each redshift. Open circles denote 
redshifts at which the VL2 catalogues are sampled while solid black lines 
trace the results derived from a cubic spline inteipolation of Rmax and Vmax- 
The dotted black lines in each panel show the expected evolution for the fit¬ 
ting functions given by|Wechsler et al.H2002) for collapse times Zc = 1-7 and 
Zc = 3.7 (see text). The VL2 curves are well contained within the shaded re¬ 
gion, which may reflect an initial “collapse” at z ~ 3.7 (in the sense of the 
Wechsler et al. formalism) followed by an episode of significant mass accre¬ 
tion at z ~ 1.7 which “resets” the concentration back to the virialization value 
of Cvir ~ 4. 


and the radius of a circular orbit of the same energy is 

GMvir 
rcirc - ■ 


(8) 


When discussing rj and e at infall we report the results of both 
methods so that we can make direct comparisons to previous 
work. 


2.5. Definition of an orbit 

In ^3.4| we compute subhalo quantities, such as tidal mass, 
taken over the course of an orbital period. To do so requires 
a precise definition of an “orbit". This is a complicated task 
since an orbit within a spherical potential is not closed, gener¬ 
ally, and traces a rosette pattern, oscillating radially between a 
mini mum pericenter, rperi. and m aximum apocenter, r^po (see, 
e.g., Binney & Tremain^|1987| l. Furthermore, the mass dis- 
tribution in realistic halos departs significantly from spherical 
symmetry, due both to triaxiality in the smooth component, 
as well as substructure. Finally, subhalos slowly spiral inward 
due to dynamical friction both from the background matter 
distribution as well as stripped material. Consequently, en¬ 
ergy and angular momentum are not in general conserved, 
and we require a robust and physical definition of an orbit that 
does not depend on simplifying assumptions such as spherical 
symmetry and conserved quantities. 

We choose to work solely from knowledge of the radial po¬ 
sition of the subhalo as a function of time, determining the 
local minima (pericenters) and maxima (apocenters). Due to 
the somewhat coarse time information, apocenters are gener¬ 
ally more accurately determined than pericenters, since ha¬ 


los spend a larger fraction of time further away from the host 
center. Thus, we dehne an orbit as that segment of the sub¬ 
halo trajectory between two successive apocenters. A given 
orbit is therefore characterized by the time at first and last 
apocenters fi and t 2 , the two apocenters rapo i and rapo, 2 . and 
the pericenter, rpeu. We take the mean of the two apocenters, 
Tapo = (rapo.i + ^apo. 2 )/ 2 , and define an effective eccentricity 


Seff = — 


while the period of the orbit is forb = f 2 ■ 


(9) 


3. RESULTS 
3.1. Host halo 

We begin by presenting the deriv ed p roperties of the host 
halo using the method outlined in §2.2| Figure shows the 
redshift evolution of the host virial mass and concentration. 
Open circles denote the redshifts for which the VL2 cata¬ 
logues are sampled while the solid black lines trace the result 
we derive after performing a cubic spline interpolation on the 
time evolution of and Vmax- Our method finds the host to 
evolve from a virial mass of Mvb =1.1 x 10" Mq at z = 4.56 
to Myii = 1.9 X 10'^ Mq at z = 0. The concentration evolves 
fr om Cvir = 2.6 at early tim es to Cvir = 12.2 at the present day. 

[Wechsler et al.| ( |2002| l showed that halo concentration is 
strongly related to mass assembly history. In particular, evolu¬ 
tion in concentration and virial mass can be described remark¬ 
ably well using a single parameter, = 1/(1 +Zc). defined as 
the formation or collapse time of the halo. They provide fit¬ 
ting relations Cvh- = 4.1a/ac and = Mvir(0)exp[-2acZ] 

which we plot in Figure [T] spanning the redshift range Zc = 
1.7-3.7. The VL2 data hts well within the shaded region 
which may reflect an initial collapse time of z 3.7 followed 
later by an episode of signihcant mass assembly at z ~ 1.7 
which resets the concentration back down to Cvir ~ 4. 


3.2. Mass functions 

The subhalo mass function provides a statistical measure 
of the amount of substructure within a host as a function of 
mass scale. In general, we can speak of two subhalo mass 
functions: the unevovled and evolved mass functions. The un¬ 
evolved mass function counts the number of subhalos based 
on their mass at the time of infall. The choice of name em¬ 
phasizes that this is a measure of the distribution of subhalos 
before they have had time to evolve under the influence of 
tidal processes within the host. The evolved mass function, 
on the other hand, counts the number of subhalos based on 
their present-day masses. 

In the middle panels of Figure]^ we plot both the unevolved 
and evolved subhalo mass functions measured from VL2. 
We compare these to the corr esponding mass functions from 
[Boylan-Kolchin et al.|(|2010| BKIO) which we re htted from 
the Millennium 11 ( Boylan-Kolchin et al.|2009 1 and Aquarius 
( jSpringel et al. 2008 1 simulations. This is a useful compari- 
son since BKlO considered host halos with similar masses to 
the VL2 host halo and use the same definition of Zinfaii as we 
do here. One difference, however, is that BKIO do not con¬ 
sider subhalos that reside outside t he vi rial radius of the host 
at z = 0 (i.e., ejected subhalos; see ^2.1| ). For more of a direct 
comparison, we also plot the VL2 mass functions with ejected 
subhalos removed. The bottom panels show more closely the 
comparison between VL2 and BKIO. 
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Fig. 2. — Top panels show the cumulative distrihution in fiQ measured at infall (left) and at z = 0 (right) for all 6145 VL2 subhalos. The gray histogram shows 
the reverse cumulative distribution function. Middle panels show the corresponding unevolved (left) and evolved (right) mass functions for all subhalos (black 
trian gles) and only “non-ejected” subhalos (gray circles). Hence, the gray circles correspond to only those 4607 subhalos currently residing within Svir at z = 0 
(see S2TJ. This is done for the purpose of compaiing to the Aquarius simulation, shown as the dark dashed blue line, based on the BKIO fitting funct ion. The 
VL2 and Aquarius mass functions show only those subhalos that survive to the present epoch. In contrast, the solid red line traces the fitting function of|Jiang &| 
|van den BoschH2014af for the unevolved mass function of all subhalos accreting onto the host. The lightly shaded curves in the left (right) panel correspond to 
evolved (unevolved) quantities in order to better show the difference between the two mass functions. The bottom left (right) panel shows the relative difference 
between the various data and the unevolved (evolved) BKIO fitting function. 


The VL2 and BKIO unevolved mass functions agree well 
with each other over most of the mass range seen here. The 
sharp cutoff at small mass simply reflects the resolution limit 
of VL2. There is considerable disagreement at the high mass 
end, though this regime is inherently noisy due to small num¬ 
ber statistics. This can be seen in the top left panel where the 
cumulative distribution in /iq at infall is shown; only 11 ob¬ 
jects with ^0 > 4 X 10“^ at infall exist. Including ejected sub¬ 
halos enhances the VL2 unevolved mass function by a con¬ 
stant factor indicating that infall mass does not play a signifi¬ 
cant role in determining whether a subhalo resides outside of 
f^vir at Z = 0. 


Note that the unevolved mass functions shown here corre¬ 
spond only to those subhalos that accrete onto the host and 
remain intac t at z = 0. The red line in Figure [2] traces equa¬ 
tion (21) of |Jiang & van den Bosch| (|2014a)l which shows 
the unevolved mass function for all subhalos ever accreted 
onto the host. This mass function is found to have a univ er¬ 
sal form (van den Bosch et ar]|2005[ [Giocoli et al.||2008] |Li| 


& Mo| 2069||, independent of host halo mass and cosmology, 


except perhaps a small dependence on tis ( Yang et al.|201 1 1. 
The main difference between this mass function and that of 
surviving subhalos occurs at high fjL where dynamical friction 
selectively disrupts massive subhalos after infall. The VL2 
unevolved mass function (black triangles) agrees well with 
the red line albeit with a small systematic shift upwards. In 
^we analyze disrupted subhalos in VL2 and find that includ¬ 
ing them here would further boost the black triangles upward 
by ~ 5% in the range 10“^ ^ Mo ^ 10“^ (see Figure 161. We 
are indeed focusing on the low-mass regime where dynamical 
friction and tidal disruption are relatively unimportant for the 
vast majority of subhalos. 

The right panels of Figure|^show the evolved counterparts. 
The evolved mass function can be thought of as a shift to 
lower mass due to tidal stripping. This can be seen by compar¬ 
ing the dark and lightly shaded blue lines. The VL2 evolved 
mass function lies systematically below the BKIO result at a 
level of about 30%. BKIO quote an intrinsic halo-to-halo scat¬ 
ter of 18% for fiQ < 10“^ which is not enough to explain the 
discrepancy seen here. Instead, the difference observed here is 
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^infall/^vir 


Fig. 3.— Cumulative (top panel) and differential (bottom panel) infall red- 
shift distributions for bins equally spaced in cosmic time. The black his¬ 
togram traces the total sample of subhalos while the blue and red histograms 
show distributions for the Icr outliers with the smallest and largest present- 
day mass ratios, respectively. 


most likely related to differences in cosmological parameters. 
In particular, VL2 uses erg = 0.74 while Aquarius simulates 
larger perturbations with ug = 0.9. It is not straightforward 
to describe how this difference manifests in the evolved mass 
function since there are at least two competing effects. On 
the one hand, the lower amplitude of fluctuations in VL2 will 
yield later formation times meaning that subhalos have less 
time on average to lose mass since infall. On the other hand, 
later formation times also yield lower subhal o conce ntrations 
which promote more efficient mass loss (see ^3.4.1 1 . 

Another factor that may contribute to this difference lies in 
the definition of tidal mass used by VL2. As described in 

E subhalo masses in VL2 are underestimated at pericenter 
; to the simplified scheme used in computing mass based 
on local density comparisons. This is in contrast to the un- 
binding procedure use d by Aquarius with the code SUBFIND 
( jSpringel et al.||^01| l. As a result, the VL2 evolved mass 
function will be biased toward smaller masses as some sub¬ 
halos will be found near pericenter at z = 0 (the unevolved 
mass function is less affected since subhalos are generally 
near apocenter at infall). Nevertheless, we only expect this to 
be a partial effect si nce a suppression of 30% was also seen in 
[Klypin et al.| ( |20ll] l when comparing VL2 and Aquarius Vmax 
functions. The physical mechanism leading to the systematic 
difference between the VL2 and Aquarius evolved mass func¬ 
tions remains to be seen. 


3.3. Subhalo properties at infall 

In this section we focus on subhalo statistics at the time of 
infall onto the host. In p articu lar, we investigate the redshift 
at which infall occurs ( |3.3.1[ ), the rad ial dis tance from the 
host at which tidal truncation initiate s ( §3.3.2| l, and show dis¬ 
tributions in orbital energy (|3.3.31l and angular momentum 
(§3.3.4 1 at infall. The results presented here are important 
as inputs into semi-analytic models of substructure evolution 


Fig. 4.— Cumulative (top panel) and differential (bottom panel) distribu¬ 
tions in the ratio of the radial distance at infall, rjnfaii, to the virial radius 
of the host, at that time. The black histogram traces the total sample 
of subhalos while the blue and red histograms show distributions for the \<j 
outliers with the smallest and largest infall redshifts, respectively. The solid 
black curve, dotted blue curve, and dashed red curve trace lognormal fits to 
the black, blue, and red histograms, re^ectively. The mean and standard 
deviation of these fits are given in Tablef^ 

and extend the results of previous works to much lower mass. 

3.3.1. Redshift: Zi„faii 

In Figure we show the cumulative distribution of infall 
redshift for dl 6145 subhalos with Zinfaii <4.56. We see that 
half of the population has fallen into the host by z = 2. We also 
plot separate distributions for the la outliers having the small¬ 
est 16% present-day mass (/i < 5 x 10“^) and largest 16% 
present-day mass (/i > 7 x 10“®). We see that presently more 
massive subhalos tend to have fallen in at more recent times. 
There are two reasons for this trend: (i) structure forms hier¬ 
archically, so halos falling in at earlier times were on average 
less massive to begin with than those infalling later; (ii) sub¬ 
halos of a given mass that fell in earlier have had more time 
to undergo tidal stripping, and will be less massive today. 

3.3.2. Radius: rinfaii 

Figurefflshows the distribution of the radial distance, rinfaii, 
between me subhalo and host at infall normalized to the virial 
radius of the host at that time. Hence, we are plotting the 
relative distance at which the subhalo has its growth history 
truncated due to tidal interactions with the host. We also plot 
separate distributions for the Icr outliers with the most recent 
infall, Zinfaii < 1-1, and earliest infall, Zinfaii > 3.93. In each 
case, the differential distribution can be well approximated 
by a lognormal form in rinfaii//?vir- The mean and standard 
deviation of the least-squared lognormal distribution for each 
population are summarized in Table 

Somewhat surprisingly, we find that over 90 per cent of sub¬ 
halos undergo tidal growth truncation outside of the virial ra¬ 
dius, with roughly 50 per cent infalling at a distance of more 
than three virial radii from the host. Considering halos falling 
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^infall 

Fig . 5.— Top panel shows the radius-redshift relationship at infall, for each 
subhalo (points), with median (solid line) and Icr distribution (shaded region) 
in each Zinf^n bin. Bottom panel shows the concentration-redshift relationship 
at infall. Solid line and and shaded region have the same meanings as in the 
top panel. Shown also is the mean concentration-redshift relationship from 
Klypinetal. (2011) for a subhalo mass of 1.4 X IO^Mq, which we find to be 
the median infall mass, independent of redshift. Halos falling in earlier are 
less concentrated and have their growth truncated at larger distances. 


in at the earliest times, Zinfaii > 3.93, this fraction rises above 
80 per cent. We emphasize, however, that these numbers are 
likely biased toward larger values since our definition of in¬ 
fall does not exclude the possibility of group preprocessing 
for which tidal truncation first occurs via interactions with ha¬ 
los other than the future host. Nevertheless, our findings are 


in qualitative agreement with past studies (Hahn et al. 


2009 


[Behroozi et al.|20T4] l showing that tidal truncation generally 
occurs outside of fJvir- This trend is also apparent in the top 
panel of Figure]^ where we show the relationship between in¬ 
fall radius and redshift directly. Why do halos at high redshift 
begin to be affected so far outside the host? 

A first hint is provided upon inspection of the bottom panel 
of Figure where the infall concentration is plotted against 
infall redshift. We see a strong correlation, with halos that 
fall in earlier having much lower concentrations. This is 
expected because the typical mass of an infalling halo does 
not change very strongly with redshift. Thus, the concentra¬ 
tions of infalling halos grow roughly as expected for halos 
of a fixed mass, e.g., 1.4 x lO^M©, which we find to be the 
median infall mass, independent of Zinfaii (see also Figures 
andp^. This is shown as the blue dashed line, which is the 
mean concentration-redshift relationship at fixed mass as de- 


TABLE1 

Lognormal ri„faii//?vir Fits. 


Population 


(T 

All subhalos 

1.09 

0.69 

Zinfall <1-1 

0.53 

0.42 

Zitifall > 3.93 

1.76 

0.61 




d 

c 



3 6 9 12 15 


Cvir 



Fig. 6.— Top panel shows the infall radius-concentration relationship at 
infall. The bottom panel shows a similar result with Cyir and rinfaii normalized 
to the concentration and scale radius of the host, respectively. In bothj)anels 
the solid line and shaded region have the same meaning as in Figurel^ 


termined by Klypin et al.| ( 201 fp 
sible explanation lies in the conci 


It would seem that a plau- 

; explanation lies in the concentration of infalling halos. 

This can be directly tested by plotting rinfaii//?vir versus Cyir, 
as shown in the top panel of Figure We see a definite 
trend of more concentrated subhalos coming closer to the host 
before undergoing infall. Our physical interpretation is that 
highly concentrated subhalos with compact density profiles 
are more resilient to tidal stripping. Hence, they are able to 
plunge deeper into the potential well of the host before ap¬ 
preciable mass loss occurs. Note that the relationship seen in 
Figure appears weaker than the trend observed when com¬ 
paring Cvir versus Zinfaii in Figure In particular, there is a 
significant fraction of halos with Cvk < 10 and r-miaw/Rvu < 3, 
with the vast majority of these falling in at late times. This 
implies there are other effects at high redshift that hinder the 
growth of infalling halos, in addition to lower central densi¬ 
ties. 

Before advancing we note that a trend in rinfaii//?vir versus 
Cvir is expected to exist even if rinfaii does not change much 
with time. This is based on the fact that both Cvir and /?vir 
generally increase with time due to the expansion of the uni¬ 


verse. To ^ to account for this, we plot, in the bottom panel 
of Figure]^ nnM\/Rs versus Cvir/Cvir, where Rs = f?vir/Cvir is 
the scale radius of the host. Normalizing Cvir and rinfaii in this 
way acts to remove redshift dependencies in Cvir and f?vir- We 
are thus plotting how close tidal truncation occurs relative to 
the central density peak of the host as a function of subhalo 
concentration relative to that of the host. Although weaker, 
we still see a definite trend of more concentrated subhalos 


® Note that the c oncentration-redshift relation observed here is systemati¬ 
cally lower than the|Klypin et al.|(2011) curve. One reason is that the latter 
was calibrated against the Bolshoi simulation which had an enhanced ampli¬ 
tude of perturbations with erg = 0.82 compared to the value of 0.74 employed 
by VL2. At a fixed redshift we expect VL2 concentrations to he lower than 
in the Bolshoi simulation based on the notion that concentration reflects the 
background density of the Universe at the time of halo formation. 
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Fig. 7.— Cumulative (top panel) and differential (bottom panel) VL2 dis¬ 
tributions in rj at the time of infall (black histogram) and first virial crossing 
with a point mass potential (gray histogram). Compared to the latter is the 
green dashed curve showing the uniform distribution reported by Z05. 

coming closer to the central depth of the host before experi¬ 
encing tidal truncation. This lends support to the notion of an 
intrinsic radius-concentration relation for infalling subhalos. 

3.3.3. Orbital energy: pmfaii 

Semi-analytic models of substructure evolution require two 
inputs as initial conditions for subhalo orbits; energy and an¬ 
gular momentum. In this section we present the infall distri¬ 
bution of energy as seen in VL2 and proceed in the next sec¬ 
tion with angular momentum. In accordance with past stud¬ 
ies, we parametrize the infall energy in terms of the variable 
77 = rdrc/Rvh, where r^rc is the radius of the circular orbit of 
the same energy as the subhalo and is the virial radius of 
the host at infall. We compute this quantitiy by first evaluating 
equation 0 for the orbital energy, and then solving equation 
0 for Tciic based on an orbit within an isolated NFW host 
potential. 

This parametrization is only valid for subhalos on bound 
orbits (E < 0). It turns out that this condition is not very re¬ 
strictive since only 38 (0.6%) of the 6145 subhalos are on 
unbound orbits at the time of infall. This number is still small 
at z = 0 when only 75 subhalos are found to be on unbound 
orbits. Interestingly, all 75 of these subhalos are outside of 
Ryir at z = 0 (i.e., they are ejected subhalos) meaning that no 
subhalos within the present virial radius are unbound. In con¬ 
trast, only a small fraction (7/38) of unbound subhalos at in¬ 
fall end up being part of the ejected population of subhalos at 
z = 0. Moreover, none of the unbound subhalos at infall are 
also unbound at z = 0. Being unbound from the host poten¬ 
tial at infall correlates neither with being presently unbound 
nor with being found outside /?vir at z = 0. Instead, it is likely 
that presently unbound subhalos acquire orbital energy from 
gravitation interactions after infall. This can be achieved, for 
example, through three-body interactions between merging 
groups of subha los as they make their first passage t ogether 
around the host ( [Sales et al.|2007[ [Ludlow et al.|2009| . 


TABLE 2 

Lognormal d//d?7 Fits. 


Population 

Zinfaii with NFW host Zeross with point host 

a cr 

All subhalos 

0.70 

0.50 

0.18 

0.40 

Z< 1 

0.38 

0.21 

0.29 

0.41 

1 <z<2 

0.53 

0.36 

0.22 

0.38 

2<z<3 

0.88 

0.44 

0.17 

0.39 

z>3 

1.50 

0.54 

0.10 

0.40 


In Figure |7] we plot the VL2 distribution in p at infall for all 
bound subhalos. This is compared to the uniform dis tribution 
between [ 0 . 6 , 1 ] used in the semi-analytic model ofjZentn^ 


et al. ( 2005 Z05) based on t he a nalysis of the N-body s imula- 
tions of|Klypin et al.|(|2001[l and [Kravtsov et al.|([2004[l. This 
also serves as the basis for t he input distribut ions of 77 used in 
the semi-analytic models of Gan et al. ( 2010[l and Jiang & van 


[den Bosch| ( [26l4b] l. The VL2 result, with a peak at 77 ~ 1.5, i 
in clear disag reem ent with the Z05 distribution. However, as 
described in |2.4| our calculation of 77 is not directly compa¬ 
rable to that of Z(j5. Firstly, Z05 report 77 at the time a subhalo 
first crosses through Ry g , whic h occurs at much later times on 
average than Zinfaii (see ^3.3.2[ l. In addition, orbital energy and 
Tcirc are computed via equations (jTj) and ([^, valid for a point 
mass host potential. The gray histogram in Figure |7] shows 
the result of applying this method to the VL2 data. The agree¬ 
ment with Z05 is better, but still heavily offset toward larger 
values of 77 . _ 

In the recent work of Jiang et al. ( 2015| J15) it was shown 
that 77 depends on mass ratio at mfall, with smaller objects 
tending toward larger 77 . A direct comparison with their re¬ 
sult is difficult, however, due to their choice of a singular 
isothermal sphere host potential. Nonetheless, their analysis 
suggests that looking only at subhalos with po 10 “^, simi¬ 
lar to those resolved in Z05, would shift the gray histogram 
in Figure [7] to the left, presumably in better agreement with 
the uniform distribution. Since VL2 samples only one host 
halo, however, we cannot test this explicitly due to insuffi¬ 
cient statistics in high-mass subhalos. Nevertheless, the result 
found here corroborates the work of J15 and indicates that 
the infall distribution assumed for 77 in various semi-analytic 
models of substructure evolution may only be strictly valid for 
relatively massive subhalos. Low-mass subhalos (p,o 10“^) 
tend to have more kinetic energy, making them less bound to 
the host, with lower specific binding energy. 

We plot in Figure [3 the dependence of 77 on mass for the 
range of mass captured in VL2. In this low-mass regime, there 
does not appear to be any significant trend of 77 with infall 
mass. Instead, we find a strong trend with infall redshift. This 
trend is attributed to the fact that J?vii is an increasing function 
of time so that 77 is pushed to larger values at earlier times. 
This trend is slightly suppressed by a competing evolution in 
Aire with redshift: we find typical values of rcirc increase by 
a factor of ~ 2 from z = 4 to z = 0 (whereas /?vir increases 
by ~ 8 in this range). This intrinsic evolution in rcirc with 
Zinfaii indicates that subhalos falling in at earlier times do so on 
slightly more bound (smaller rdre) orbits. The top right panel 
of Figure [^ shows a strong trend in 77 with Cyir which can be 
attributed to the concentration-redshift relation seen earlier in 
Figure]^ 

We note that the trends observed in the top panels of Figure 
[^depend on which version of 77 we present. Namely, we find 
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Fig. 8.— The dependence of rj (top panels) and e (bottom panels) at infall on mass ratio at infall (left), infall redshift (middle), and concentration at infall 
(right). In both cases we assume an NFW potential for the host. Points denote individual subhalos while the solid line shows the median trend for bins along the 
jc-axis with the corresponding Icr spread indicated by the shaded region. 


that using r] computed at virial crossing for a point mass host 
potential mostly washes out the dependence of rj with redshift 
and concentration. The reason is that subhalos tend to become 
more bound in time after infall so that rj is pushed to smaller 
values at the time of virial crossing, Zcross. which generally 
occurs after Zinfaii (see Figure]^. This seems to occur in such 
a way as to mostly cancel the evolution in /?vir with z- The 
difference between the black and gray histograms in Figure 
can therefore be mostly explained by two effects; (1) larger 
f?vir Across and (2) Smaller rat Zcross- 
We find that the probability distributions, d//dry, can be 
well modelled by lognormal distributions. In Table we list 
the best-fit mean and standard deviations for both me total 
distribution as well as those obtained from different redshift 
ranges. The latter corresponds to cuts in Zinfaii and Zcross for 
the infall and virial crossing methods, respectively. These fit¬ 
ting functions are appropriate for subhalos with mass ratio 
Mo ^ 10“^. Higher mass subhalos should shift closer to the 
Z05 curve in Figure |7] In ^we find that subhalos on tightly 
bound orbits with r] < I are preferentially disrupted after in¬ 
fall. However, this has only a small impact on the distributions 
presented here since this bias is small and there are far fewer 
disrupted than surviving subhalos. Hence, the fitting func¬ 
tions provided here should be applicable to the total ensemble 
of subhalos (surviving plus disrupted) that ever accreted onto 
the host. 


3.3.4. Orbital angular momentum: ei„faii 


Stu dies of substructure evolution (e.g., Penarrubia et al. 
|2008| l show that subhalos on more radial orbits with lower 
specific angular momentum plunge deeper into their hosts and 
experience accelerated mass loss over subhalos on more cir- 
cular o rbits with higher specific angular momentum (see also 
^3.4.1|). Accurately modelling subhalo evolution therefore re¬ 
quires a good handle on the distribution of angular momentum 
at the time of infall. As such, a great deal of work has been 


done on measuring this distribution from N-body simulations 


(Navarro et al.|1995 Tormen 

l9971|Ghigna et al.|l998 |Ben- 

son|2005HZentner et al.|2005 

Wang et al.|2005[|khochfar & 

Burkert 

sions of 

|2UU6 Jiang et al.||2UU8 Wetzel|j2Ul 1||. The conclu- 

these works agree well with each other; the circular- 


ity distribution of infalling satellites is peaked at e « 0.5 and 
falls off on either side so that neither largely radial (e ~ 0) 
nor largely circular (e ~ 1) orbits occur. Below we turn our 
attention to the circularity distribution measured in VL2. 

In Figure]^ we plot the infall distribution in e for all bound 
subhalos. As in Figure [7] we show the result at Zinfaii for an 
isolated NFW host potential (black histogram) as well as the 
result at first /?vir crossing for a point mass host potential. The 
latter can be compared to the various curves showing the in¬ 
fall distributions used in semi-analytic models of substructure 
evolution. First , the blue dotted l ine is the Gaussian distri¬ 
bution used by Taylor & Babul ( |2004| l with mean e = 0.4 
and standard deviation cr = 0.28 which was selecte d so that 
the fin al di stribution at z = 0 ma tches the results of ITormenl 
(|1997|) and|Ghigna et al.| (|1998||. Second, the green dashed 
line shows the one-parameter fi distribu t ion use d in the mod¬ 
els of Z05 and |Jiang & van den Bosch| (|2014b|). Finally, the 
red dot-dashed curve sho ws the infall distrib ution assumed in 
the semi-analytic model of |Gan et'ar] ( |2010] l, which was taken 
from the analysis of the hydrodynamic simulations of IJiangl 
|etar| ( |2008l ). 


The gray histogram, with median e = 0.55, is in reasonable 
agreement with the various curves used in semi-analytic mod¬ 
els. There is, however, a clear excess in nearly circular orbits 
with e ~ 1. The reason for this is the same as was discussed 
in the previous section. It was shown in JI5 that circularity 
is highly dependent on mass, with high mass subhalos tend¬ 
ing to move along radial orbits while low mass subhalos tend 
to have more circular orbits. One possible explanation is re¬ 
lated to the environment in which these objects form. High 
mass halos are more biased towards forming in high density 
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Fig. 9.— Cumulative (top panel) and differential (bottom panel) distribu¬ 
tions in e at the time of infall (black histogram) and first virial crossing with 
a point mass potential (gray histogram). Compared to the latter are various 
curves showing the circularity distributions employed in semi-analytic mod- 
els of substructure evo lution: the blue dotted curve is the distribution from 
ITaylor & Babul|j2004f; the green dashed curve is the distribution given by 
Z,U5; the red dot-dashed curve is the distribution from|Gan et al.|(2010). 


regions such as filaments and are consequently more likely 
to fall radially into their host with lower specific angular mo¬ 
mentum. Low mass subhalos are less biased to forming within 
filaments and are thus more likely to fall into their host with 
a larger component of tangential motion. Another possible 
explanation, independent of the first, is that low-mass subha¬ 
los are simply more likely to acquire tangential motion from 
gravitational interactions with nearby massive objects prior to 
infall. 

In the bottom panels of Figure we show the dependence 
of e on mass, redshift, and concentration. As with rj, we find 
that e does not exhibit strong mass dependence in the range 
Fo ^ 10“^ probed by VL2. There also does not appear to be 
m uch dependenc e on Zinfaii, consistent with the earlier work 
of |Wetzel1 ( |201 l| l. There is correspondingly little dependence 
on e with infall concentration. Replacing the lower panels of 
Figure]^ with e computed at virial crossing for a point mass 
host potential produces the same (lack of) trends. 

Even though the two methods produce similar trends in Fig¬ 
ure it remains difficult to qualitatively describe the differ¬ 
ences between the black and gray histograms in Figure]^ The 
reason is the dependence of on the form of the host poten¬ 
tial assumed. In fact, we find that e, unlike rj, is very sensitive 
to the form of the host potential. For example, computing e 
at infall with a point mass host potential pushes the black his¬ 
togram in Figure]^ to even larger values. Counter intuitively, 
computing e at virial crossing with an NFW host potential 
also pushes the gray histogram to larger values, being almost 
on top of the black histogram. We therefore note that one 
should be careful in choosing a fitting function for e that suits 
their specific needs. 

We fit the circularity distributions in Figure]^ with the fol¬ 


Fig. 10.— Cumulative (top panel) and differential (bottom panel) distri¬ 
butions in the mass loss over the course of an orbit normalized to the mass 
at infall. The black histogram pertains to the first orbit after infall while the 
blue histogram shows the last orbit before z = 0. Only those 2714 subhalos 
(44% of the total population) that complete at least two orbits after infall ai‘e 
shown. 


lowing form: 

'^=ae°‘(b-ef. (10) 

de 

We find the best-fit coefficients (a, b, a, /3) = (3.696, 1.12, 
1.07, 0.68) at infall with an NFW host potential and (a, b, a, 
P) = (1.508, 1.77, 1.05, 2.45) at virial crossing with a point 
mass host potential. These fitting functions are appropriate 
for subhalos with mass ratio /tq ^ 10 “^ and are independent 
of redshift. Higher mass subhalos should have distributions 
in closer agreement with the other fitting functions plotted in 
Figure!^ In ®we show that circularity does not influence the 
surviv^ility of low-mass subhalos. The fitting functions pro¬ 
vided here are thus applicable to the total ensemble of subha¬ 
los (surviving plus disrupted) that ever accreted onto the host. 


3.4. Evolution 

In this section we focus on the evolution of subhalo proper¬ 
ties over the course of infall to the present day. This includes 
internal subhalo properties such as tidal mass and central den¬ 
sity as well as orbital properties including radial period and 
angular momentum. Our results are used to test some of the 
fundamental assumptions underlying models of substructure 
evolution. 

A common assumption in modeling tidal mass stripping 
is that subhalos of a given mass lose a certain fr action of 
their mass in one dynamical time. For example, van den 
|Bosch et al.| ( |2005l l developed a model in which the mass loss 
rate of a given subhalo is m cx where r(iyn(z) oc 

(Apciit(z))~^^^ is proportional to the free fall time of a halo, in- 
dependent of mass. Recently, |Jiang & van den Bosch|(|2014b i 
used numerical simulations to ht the mass dependence, hnd- 
ing C, = 0.07. This is very close to the case C = 0, in which the 
fractional mass loss rate is independent of mass. 

































Evolution of Low Mass Galactic Subhalos and Dependence on Concentration 


11 





Fig. 11.— Left (right) panels show the fractional amount of mass, Am, lost over the course of the first apocenter-to-apocenter orbit after infall (last apocenter- 
to-apocenter orbit before z = 0) normalized to the mass, mi„faii, of the subhalo at the time of infall. The various rows show: (i) fractional mass lost as a function 
of m/Mvir with both quantities computed at the start of the orbit; (ii) fractional mass lost as a function of Cmax/Cmax with both quantities computed at the start of 
the orbit; (iii) fractional mass lost as a function of eccentricity Cf-ft computed from equation |^. In each panel the solid black line traces the median mass loss for 
bins with equal number of samples in the x axis while the shaded region shows the la spread about this line. In the top and bottom rows, blue circles and red 
triangles highlight the 2a outliers with the smallest 2.3% and largest 2.3% values of Cmax/Cmax. In the middle row, blue circles highlight the 2a outliers with the 
smallest eccentricities while red triangles highlight the 2a outliers with the largest eccentricities. 


Modeling dynamical friction and sinking due to the resul¬ 
tant loss of angular momentum plays a prominent role in mod¬ 
eling subhalo orbital evolution. The trajectory of subhalos 
through the host must be modeled accurately. Assumptions 
typically involve spherical symmetry, wherein the torque of 
dynamical friction is in the direction of the subhalo orbital 
angular momentum, and subhalos orbit in the same plane. 

The question naturally arises whether common assumptions 
such as those discussed above hold for the low-mass subhalos 
considered here. In the following sections, we examine sepa¬ 
rately the mass lost per orbit per halo mass, the orbital period 


per host dynamical time, and the alignment of tidal torques 
and angular momentum. 

3 . 4 . 1 . Tidal mass loss: Am/m 

The general picture of mass loss is related to the processes 
of dynamical friction and tidal stripping, as follows. The con¬ 
tinued force of dynamical friction causes an infalling subhalo 
to slowly descend into its host. As the orbital radius shrinks, 
so too does the tidal radius, causing the subhalo to continu¬ 
ally shed mass from the outside-in. The internal structure of 
the subhalo is also affected, generally puffing outwards due 
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to the injection of tidal heat, promoting additional mass loss. 
Mass loss will vary over the course of an orbital period, being 
strongest (weakest) at pericenter (apocenter) when tidal inter¬ 
actions with the host are greatest (smallest). Tidal mass loss 
is clearly a complicated process that will depend on both the 
internal structure of a subhalo as well as its orbital parameters. 

We begin our investigation of mass loss in Figure [T0|where 
we plot distributions in Am/minfaii. Here, the mass change 


Am = , 


i,i“'«apo .2 where m^po.i and mapo ,2 are the mass 


at the start and end of the orbit, respectively. The black his¬ 
togram shows mass loss over the course of the first orbit after 
infall while the blue histogram shows mass loss over the last 
orbit before z = 0. Recall that we define an orbit to corre¬ 
spond to the time between successive apocenter passages so 
the “first” orbit does not start exactly at infall and the “last” 
orbit does not end exactly at z = 0. More specifically, the first 
orbit begins at the first apocenter after infaljj while the last 
orbit terminates at the last apocenter before the present time. 
We find that 3966 (65%) subhalos finish at least one orbit af¬ 
ter infall while 2714 (44%) finish at least two. Since we are 
interested in comparing how mass loss changes with time, we 
plot only those 2714 subhalos for which the first and last orbit 
is different. 

Comparing the two distributions in Figure [TOl shows that 
subhalos tend to lose a larger fraction of their initial mass dur¬ 
ing their first orbit compared to their last orbit. In particular, 
the median mass loss in the first orbit is 16% of the initial 
mass while the median mass loss in the last orbit is about an 
order of magnitude smaller, at 3% of the initial mass. Note 
that not all subhalos lose mass over the course of an orbital 
period. In particular, for both the first and last orbit, roughly 
5% of subhalos actually gain mass. This likely occurs either 
through direct merger with smaller systems or, more gradu¬ 
ally, through the accretion of surrounding material. 

We proceed to investigate the dependence of mass loss on 
subhalo properties. The top row of Figure!^ shows mass loss 
versus mass ratio at the start of the orbit. FOTboth orbits, more 
massive subhalos tend to lose more mass on average. Nor¬ 
mally, we would expect this result on the basis of a dynamical 
friction argument whereby the oribts of massive subhalos are 
preferentially dragged into the depths of the host, promoting 
enhanced mass loss. However, we do not expect this argument 
to apply here since the dynamical friction merging timescale 
for jj, < 10~^ subhalos is mu ch longer than the Hubble time 
(|Boylm-Kolchin et al.|2008|l. 

Instead, the observed correlation with mass is the result of 
the mass-concentration relation which states that more mas¬ 
sive subhalos will be less concentrated on average. We plot 
as blue circles (red triangles) the 2tT outliers with the small¬ 
est 2.3% (largest 2.3%) values of Cmax/Cmax at the start of the 
orbit. From the definition of Cmax in equation Q, the ratio 
Cmax/Cmax describes the relative central density ofthe subhalo 
to the host. In both panels a clear dichotomy emerges with the 
least (most) concentrated, and most (least) massive, subhalos 
loosing (retaining) more mass per obit. 

This is made more apparent in the middle row of Figure 


11 where we see a strong negative slope in mass loss versus 


concentration. There is still considerable scatter at fixed con- 


^ We find that infall is roughly symmetric about the turnaround point where 
the subhalo first detaches from the Hubble flow and begins its descent towards 
the host. In particular, 52% of subhalos start to lose mass before turnaround 
while 48% begin losing mass after turnaround. Hence, for roughly half of the 
cases, the first orbit begins at the turnaround radius, corresponding to the first 
apocentre. 



Fig. 12.— Orbital period, rorj,, in units of the halo dynamical time, Tjyn, for 
subhalos on their first orbit, colored by concentration, as a function of first 
apocenter, rapo, in units of the host virial radius. All quantities are determined 
at the beginning of the orbit. Shown as the dashed line is the relationship 
expected for radial orbits (eccentricity e= 1) orbiting a point with mass Myip. 


centration which can be partly attributed to eccentricity. Com¬ 
paring blue circles and red triangles shows that for fixed con¬ 
centration, more circular (radial) orbits tend to retain (lose) 
more mass on average. A direct comparison is plotted in the 
bottom row of Figure where we detect a small correlation 
between mass loss and^centricity. 

We have also checked for correlation between mass loss and 
pericenter, rpeu. One would expect that subhalos plunging fur¬ 
ther into the depths of the host, where tidal forces are strong, 
would experience enhanced mass loss. We instead find almost 
no correlation with r-perj^ The reason is that subhalos closer 
to the host center tend to be more concentrated (see g3.3.2| i 
which washes out the dependence on rpeu. 

We conclude that tidal mass loss in the regime of low-mass 
subhalos is most directly correlated with concentration. When 
concentration is held fixed, we find no trend in mass loss with 
varying mass. The apparent trend seen when comparing mass 
loss versus mass is simply a reflection of the fact that mass 
is correlated with concentration. This result makes physical 
sense in the limit of weak dynamical friction since it is the 
density of a subhalo, relative to its host, that determines how 
tightly a subhalo o n a stable orbit retains its contents (e.g., 
[Taffoni et al.|2003 1. At fixed concentration, subhalos on more 
eccentric (i.e., radial) orbits tend to lose more mass than sub¬ 
halos on circular orbits. This may highlight the importance 
of tidal heating which results when a rapidly varying gravi- 


^ Note that our calculation of Tperi is based upon a cubic spline interpo¬ 
lation of subhalo radial distance from discrete VL2 snapshots (see The 
snapshots are separated by 0.688 Gyr which may lead to a potentially crude 
estimation of the true rped- Another way to estimate rperi is to solve the roots 
in the equation of motion of a point particle of e nergy E and angular momen - 
tum L in a static NFW host potential (see, e.g., |Binney & Tremaine|1987^ . 
This may lead to a better determination of rperi E and L vary more snioothly 
with time than radial distance. Nevertheless, we have checked that using rperi 
computed in this way changes neither the results on mass loss versus eccen¬ 
tricity nor mass loss versus pericentre. 
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Fig. 13.— Cumulative (top panel) and differential (bottom panel) distribu¬ 
tions in the dot product between the angular momentum normal vector at z = 0 
and infall. The black histogram traces the total sample of subhalos while the 
blue and red histograms show the Icj outliers with latest and earliest infall, 
respectively. Those subhalos that fell in recently tend to have their present- 
day angular momentum vector more closely aligned with the infall direction 
while subhalos that fell in earlier approach a uniform distribution in x (shown 
for comparison as the dotted red line). 

tational potential injects energy into subhalos, 
outwards and promoting further mass loss (e.g., 

[ 200 ^ . 

3.4.2. Orbital period 

The dynamical time for a halo is usually defined as the free- 
fall time of a test particle in a static, uniform sphere at the 
virial density, 

r2^„^[16GperitA/(37r)r' = ^^. (11) 

It is natural to expect, all else being equal, that timescales 
within the halo should scale in proportion to this dynamical 
time. For example, the time to complete an orbit is roughly 
proportional to the dynamical time. Similar scaling arguments 
apply for timescales other than orbital period, such as the tidal 
mass loss time, mjrh. Departures from a simple linear scaling 
with the dynamical time occur because orbital shapes vary 
from subhalo to subhalo. For example, halos on larger orbits 
should have longer orbital times, with a correlation between 
semi-major axis and period that reflects the mass distribution 
of the host halo around the virial radius. 

In Figure we show the orbital period of halos just after 
they fall in, defined as the time between the first two apocen- 
ters after infall. The orbital period is plotted in units of the dy¬ 
namical time, Tdyn, at the beginning of the orbit. On the x-axis 
we show Tapo/f^vir- The first striking trend is the large spread 
in forb/Tiyn values. There is also a similarly large spread, about 
an order of magnitude, in the apocenters, rapo/f^vir of halos on 
their first orbit. In fact, there is a strong correlation between 
orbital period and apocentric distance, as expected. Interest¬ 
ingly, there are two “clouds" of subhalos. Those with high 


concentration at small radius, and those with low concentra¬ 
tion at large radius. 

Also shown is the relationship expected for radial orbits 
(where the semi-major axis a = rapo/2) around a point mass 
with M = Mvir, = 47r2aV(GA/vir) = 7 rV3p^/(2GMvir). Com¬ 
bining with equation ( [TT] l, we obtain 

(12) 

Ttyn \ 2/Vvij. J 

Subhalos would lie along this line only if they were on radial 
orbits and all the halo mass was located at the halo center. 
In general, departures from radial orbits {a > r^^o/2) lead to 
longer orbital times, while the presence of matter outside the 
virial radius (i.e. the overdensity associated with continuous 
infall onto the host) leads to shorter orbital times. The latter 
effect could be responsible for the shorter times at rapo > 
although more information about the evolving density profile 
outside the virial radius would be required to make a quanti¬ 
tative comparison. 

3.4.3. Orbital plane 

A common assumption made in models of substructure evo¬ 
lution is spherical symmetry of the host. Subhalo orbits are 
generally integrated in either a static potential or one that dy¬ 
namically adjusts (e.g., through mass accretion) in a spheri¬ 
cally symmetric manner. In either case, the direction of the or¬ 
bital angular momentum vector is conserved since it is aligned 
with the direction of the torque. Hence, an obvious test of 
spherical symmetry within VL2 is to look for changes in the 
orientation of the orbital plane. 

In Figure [T^ we plot the distribution in the dot product be¬ 
tween the angular momentum normal vector at z = 0 and infall; 

X “ To ■ Tinfall. (13) 

The black histogram shows the distribution for all subhalos 
while the blue and red histograms show the Icr outliers with 
the latest and earliest infall, respectively. The median x for 
all subhalos is 0.39 while subhalos with the earliest and latest 
infall time have median values 0.16 and 0.88, respectively. 
There is a clear trend of recently infalling halos remaining in 
the same orbital plane while subhalos with early infall have 
their orientation randomly aligned. 

Subhalo orbits are continuously torqued after infall, in a di¬ 
rection that is not aligned with the angular momentum vector. 
Subhalos spending more time in the host experience larger 
changes in L. 98% of subhalos with Zinfaii < 11 do not fin¬ 
ish an orbital period by z = 0 while 65% of subhalos with 
Zinfaii > 3.93 finish at least three orbits. The latter population 
approach a uniform distribution in x, indicating that memory 
of the initial orbital plane is lost after a few orbits within the 
host. It is clear that the assumption of spherical symmetry 
does not apply. 

This result is not too surprising, however, since dark matter 
halos are generally triaxial in shape and the host will expe¬ 
rience anisotropic mass redistribution as massive o bjects are 
bi ased to ward filamentary accretion (see also, e.g., |Zemp etj 
|ar||2009| ). Another possible source of or bital torque is sub- 
structure interaction. [Slater & Bell P013] l used VL2 to show 
that a significant fraction of subhalos accrete as groups with 
correlated trajectories that lead to frequent interaction over 
time. Such interactions can lead to a complex redistribution 
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Fig. 14.— Evolution of Vmax (top), rmax (middle), and pmax (bottom) as a 
function of mass retained at ^ = 0. In each panel, points represent individual 
subhalos coloured according to concentration at infall while the solid black 
line shows the median trend with the associated Icr scatter shaded in gray. 
The PIO relation for isolated NFW subhalos of fixed concentration falling 
into a static host is shown as the dashed black line in each panel. The thin 
da shed gray line in the top panel tr aces the Bolshoi Vmax relation as reported 
injvan den Bosch & Jiang] 42014f while the dashed purple line traces the 
corrected relation including only sufficiently resolved subhalos in Bolshoi 
(van den Bosch, private communication). 

of orbital energy and angular momenta for the l ow-mass sub¬ 
halos that are abundant here ( [Sales et al.|2007t [Ludlow et al.j 
|2009| l. Though a more detailed inspection of orbits is required 
to assess the significance of these effects in VL2, our result 
highlights the importance of considering host anisotropy and 
subhalo interaction in semi-analytic models of substructure 
evolution. 

3.4.4. Subhalo internal structure: rmax and Vmax 

The main observable properties of luminous subhalos are 
their velocity structure, often described in terms of the cir¬ 
cular velocity profile, = GM(< r)/r. In particular, most 
dynamical measurements provide robust constraints on the 
maximum circular velocity, Vmax, and the radius at which this 
occurs Tmax- In this section we show evolution in these two 
quantities as subhalos descend into the host. We refer to this 
as evolution in internal structure in the sense that Vmax and 
r„,ax describe central density with p„,ax oc (Vn,ax/?‘max)^ being 
the mean density within r^ax. _ 

A number of previous works (IHayashi et al.[2003[[Penarru- 
[bia et al.|2008| 2010[ ) have studied the evolution in Vmax and 
Vmax using numerical simulations where isolated subhalos are 


dropped into the potential of a static host. These studies come 
to the same conclusion that Vmax and rmax evolve along tightly 
defined trajectories when written in terms of the mass fraction 
retained after infall. In particular, defining x= m/minMi and 
taking y to represent either rniax/Wanfaii or Vmax/vmaxinfaii, it 
is found that subhalos starting at (.r,y) = (1,1) move steadily 
along the track 

where a and /3 are fitting coeffic ients. For the case of NFW 
subhalos, [Penarrubia et al. ( 2010[ PIO) find the result (a,f3) = 
(0.4,0.3) for Vmax and (a,P) = (-0.3,0.4) for rmax- 

In Figure [T4| we plot the ratio of the present-day values of 
Vmax. ^max. and pmax to their infall values versus the fraction 
of mass retained at z = 0. In each panel, the solid black line 
traces the median relation and the shaded region shows the Icr 
scatter. The vast majority of subhalos experience reduction in 
^max and Vmax. with a larger suppression in the former, lead¬ 
ing to a net increase in pmax with increasing mass loss. The 
median relation in each panel can be compared to the dashed 
black line showing the PIO result. We find VL2 agrees well 
with PIO for Vmax but begins to diverge at low mass reten¬ 
tion for r^iax and Pmax- The dashed gray line in the top panel 
shows eq uation ( 14 1 with (a,/?) = (0.60,0 .44) which was re¬ 
ported by [van den Bosch & Jiang[ ([2014j| to ht evolution in 
Vmax for subhalos in the Bolshoi simulation. The VL2 data 
sits systematically above the Bolshoi result. The dashed pur¬ 
ple line shows a corrected form (a, /3) = (0.36,0.33) which fits 
the Bolshoi relation when insufficiently resolved subhalos are 
removed from the sample (van den Bosch, private commu¬ 
nication). This shows much better agreement with the VL2 
result. 

The points in Figure are coloured in terms of subhalo 
concentration at infall. In the case of Vmax we do not see much 
dependence on concentration other than the fact that subha¬ 
los with larger Cvir tend to have fallen in more recently and 
therefore have not had as much time to evolve to the left side 
of the plot. In contrast, rmax and pmax show strong stratifi¬ 
cation in Cvir with the least concentrated subhalos showing 
systematically greater reduction in rmax and enhancement in 
Pmax- We offer a heuristic explanation as follows. Subhalos 
on slowly sinking orbits experience mass loss until the tidal 
radius shrinks to the point at which the mean interior density 
is proportional to the local density of the host. Since subha¬ 
los are exposed to (roughly) the same local density, those that 
were initially more dense (i.e., larger Cvh) naturally approach 
a smaller value of Pmax/Pmax.mfaii at late times. 

PIO do not find significant scatter since they consider sub¬ 
halos of fixed concentration. They do find, however, that vary¬ 
ing the shape of the subhalo inner density profiles at fixed con¬ 
centration changes the coefficients a and /3. We find the com¬ 
plementary result that varying concentration at fixed shape 
leads to substantially different structural evolution. 

This point is made more illuminating by plotting Vmax ver¬ 
sus Tmax. as in Figure In the left panel we colour points 
according to the mass retained at z = 0 while the right panel 
shows concentration at infall. The left panel shows a clear 
gradient in colour, reinforcing the notion of prev ious works 
(Hayashi et al. [2003[[Penarrubia et al.[2008||2010|l that evolu- 
tion in internal structure does not depend on how mass is lost, 
but only how much mass is lost. As expected, however, subha¬ 
los are not bound to a single trajectory in the Vmax“^max plane. 
The right panel shows that scatter at fixed mass loss can be 





































Evolution of Low Mass Galactic Subhalos and Dependence on Concentration 


15 




Fig . 15.— Evolution of subhalos in the normalized Vmax - ^max plane. Left panel shows subhalos coloured according to the fraction of mass remaining at z = 0 
while the right panel colours subhalos accordi ng t o their concentration at infall. The thick black line shows the median trend with the associated Icj scatter shaded 
in gray. The thin purple line shows equation i |14| with (q;,/ 3) = (0.25,0.34) chosen to match the median relation. The dashed black line shows the PIO relation 
for isolated NFW subhalos of fixed concentration falling into a static host. 


attributed to concentration, as in Figure [T4| The median evo¬ 
lution for all subhalos is shown as the solid black line in each 
panel with Icr scatter shaded in gray. The purple line fits the 
median trend using equation ( [l4| ) with (a,/?) = (0.25,0.34). 
This sits above the PIO result that was derived from subhalos 
of fixed Cvir = 23; a considerably larger value than the median 
concentration of 7 found in VL2. 

4. DISRUPTED SUBHALO POPULATION 

Up to this point we have only considered the population 
of subhalos that survive to the present day. This leaves open 
questions regarding any potential biases that may exist in our 
results due to the omission of disrupted subhalos. In this sec¬ 
tion we attempt to address these questions by making use of 
the second VL2 public catalogue. 

The second catalogue contains a similar set of evolution¬ 
ary tracks as the main catalogue except that it pertains to 
the 20000 largest systems in the simulation box identified at 
z = 4.56. The two catalogues are not mutually exclusive as 
some of the surviving subhalos at z = 0 also happened to be 
of the largest systems present at z = 4.56. The utility of the 
second catalogue is that it contains subhalos that disrupt prior 
to the present day. Though this does not constitute the full 
ensemble of disrupted subhalos (some subhalos that disrupt 
were not of the largest systems at z = 4.56) it should be enough 
to elucidate differences between surviving and disrupted sub¬ 
halos. For convenience we henceforth refer to the main VL2 
data set (considered in all preceding sections) as VL2-D1 and 
refer to the second data set as VL2-D2. 

We apply the same framework outlined in ^to VL2-D2. 
Namely, we identify subhalos as those systems that at some 
point passed within the instantaneous virial radius of the host. 
Of these subhalos we identify the surviving population as 
those that still exist as intact objects at z = 0. Conversely, the 
disrupted population consists of those subhalos that fall below 
the mass resolution of the VL2 halo finder some time before 
reaching z = 0. We find a total of 4843 subhalos in VL2-D2 
of which 3843 (79%) belong to the surviving group and 1000 
(21%) belong to the disrupted group. Note that ~ 63% of the 


6145 subhalos from VL2-D1 are also part of VL2-D2. 

We examine in Figure[T^the dependence of survivability on 
infall redshift, mass ratio, and concentration. In the top pan¬ 
els, we plot the distribution in each quantity for all surviving 
(disrupted) subhalos from VL2-D2 as solid blue (dashed red) 
lines. For comparison, the solid black line traces the distribu¬ 
tion in each quantity for the VL2-D1 subhalos. The vertical 
dotted lines denote the median value of the distribution with 
the corresponding colour. In the bottom panels, the dashed red 
line traces the ratio, Nd 2 (D)/Nu 2 {S), of disrupted to surviving 
subhalos in VL2-D2, with the horizontal red line denoting the 
mean value 1000/3843. To test for bias in the surviving frac¬ 
tion of VL2-D2 subhalos, the solid blue line traces the ratio, 
Nd 2 (S)/Ndi, of the surviving subhalos in VL2-D2 to the full 
ensemble of surviving subhalos contained in VL2-D1. The 
horizontal blue line denotes the mean fraction 3843/6145. 

We begin with infall redshift. The blue line in the bottom 
left panel shows that VL2-D2 is slightly biased towards con¬ 
taining those surviving subhalos with larger values of Zinfaii- 
This reflects the fact that subhalos with smaller values of Zinfaii 
were less likely to exist as large objects at z = 4.56 when VL2- 
D2 was constructed. Comparing the surviving and disrupted 
populations in VL2-D2 shows that subhalos with Zinfaii ^ 3 are 
much more likely to belong to the latter group. That is, sub¬ 
halos spending more time exposed to the tidal field of the host 
are more likely to disrupt by the present day. 

Next we investigate mass ratio at infall. The middle pan¬ 
els of Figure [T^show little dependence of survivability in the 
mass range 10“^ ^ ^ 10“^. For larger masses, we expect 

the disrupted fraction to increase as dynamical friction prefer¬ 
entially causes massive subhalos to plummet into the depths 
of the host where tidal forces are strongest. We indeed see an 
upturn at /To ^ 10“^ though VL2 is hindered by small num¬ 
ber statistics in this regime to make a meaningful statement 
here. We also observe a rapid rise in the disrupted fraction for 
do ^ 10“^. This is expected due to the finite mass resolution 
of the simulation - subhalos closer to the resolution limit are 
more likely to “disrupt”. There is a related drop in the blue 
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Fig. 16.— Top panels show the distribution of subhalos based on Zinfeii (left), mass ratio at infall (middle), and concentration at infall (right). The solid black 
line shows the distribution of subhalos from the main catalogue, VL2-D1, considered in the preceding sections of this paper. The solid blue (dashed red) line 
traces the distribution of surviving (disrupted) subhalos from VL2-D2. The vertical dotted lines show the median value of the distribution of the corresponding 
colour. The dashed red lines in the bottom panels trace the ratio A^d 2 (D) /M) 2 (S) of disrupted to surviving subhalos in VL2-D2 with the dotted red line showing the 
mean value 1000/3843. The solid blue lines trace the fraction Wd 2 (S)/^di of surviving subhalos in VL2-D2 to surviving subhalos in VL2-D1. The dotted blue 
line shows the mean fraction 3843/6145. 


line in the bottom middle panel indicating that the missing 
surviving population in VL2-D2 are almost exclusively low- 
mass subhalos. This reflects the early infall bias in VL2-D2: 
subhalos nearer the resolution limit are only likely to survive 
to z = 0 if they infall later. 

We finish by examining infall concentration in the right 
panels of Figure [T^ First, we see that VL2-D2 is biased to¬ 
wards containing low-concentration surviving subhalos, con¬ 
sistent with the early infall bias via the concentration-redshift 
relation (see Figure 1^. Moreover, of all subhalos in VL2- 
D2, those with Cvir ^5 are much more likely to disrupt by 
z = 0. This is consistent with our previous finding that low- 
concentration subhalos are more susceptible to tidal stripping 
from the host (see Figure [TT]). In fact, the strong dependence 
of survivability on Cvir and not on iiq strengthens the notion 
that mass loss is more strongly connected to concentration 
than mass for low-mass (po ^ 10 “^) subhalos. 

We now shift focus to the infall distribution of orbital en¬ 
ergy and angular momentum for surviving versus disrupted 
subhalos. In Figure [T7| we plot the infall distribution of rj for 
the surviving (disrupted) subhalos in VL2-D2 as a blue (red) 
histogram. The blue histogram can be compared to the black 
histogram showing the infall distribution for all surviving sub¬ 
halos in VL2-D1 (i.e., the black histogram in Figure!^. We 
see good agreement between the surviving VL2-D2 sShalos 
and the full VL2-D1 population. The disrupted population, on 
the other hand, agrees well with the other curves for p > 2 , but 
displays a much flatter distribution for smaller p. Most impor¬ 
tantly is the excess at 77 < 1 , indicating that subhalos strongly 
bound to the host at infall are preferentially disrupted by z = 0 . 

This result suggests that the black histogram is suppressed 
at p < 1 compared to the distribution of all subhalos that ever 
fell onto the host, regardless of survivability. The gray his¬ 
togram in Figure shows the result of combining the dis- 



Fig. 17.— Cumulative (top panel) and differential (bottom panel) distri¬ 
butions in 77 for all surviving VL2 subhalos (black histogram), the surviving 
subhalos in the VL2-D2 catalogue (blue histogram), and the disrupted subha¬ 
los in the VL2-D2 catalogue (red histogram). The gray histogram shows the 
result of combining the total population of surviving subhalos with the VL2- 
D2 disnapted population. In each case we report r] at Zinfaii for NFW host 
potential. The results are qualitatively similar if we consider virial crossing 
with a point mass host potential. 

rupted subhalos in VL2-D2 with the VL2-D1 catalogue. Do¬ 
ing so results in only a minor change to the black histogram 
since there are six times fewer disrupted subhalos than sur- 
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Fig. 18.— Cumulative (top panel) and differential (bottom panel) distri¬ 
butions in e for all surviving VL2 subhalos (black histogram), the surviving 
subhalos in the VL2-D2 catalogue (blue histogram), and the disrupted sub¬ 
halos in the VL2-D2 catalogue (red histogram). In each case we report e at 
Zinfaii for an NFW host potential. The results are qualitatively similar if we 
consider virial crossing with a point mass host potential. 


viving subhalos. Though we do not have access to the full 
ensemble of disrupted subhalos, it seems unlikely that we are 
missing a large enough fraction for there to be a significant 
impact on the inf all distribution and fitting functions for 77 pre¬ 
sented in §3.3.31 

In Figure 18 we plot the infall distribution of e for the sur¬ 
viving (disrupted) subhalos in VL2-D2 as a blue (red) his¬ 
togram. The surviving population can be compared to the 
black histogram showing the infall distribution for all surviv¬ 
ing subhalos in VL2-D1 (i.e., the black histogram in Figure]^. 
In this case, we see good agreement between all histograms. 
Hence, unlike orbital energy, angular momentum does not ap¬ 
pear to play a significant role in determining the survivability 
of subhalos. T he infa ll distribution and fitting functions for 
e presented in §3.3.4| are therefore robust to the inclusion of 
disrupted subhalos! 


5. SUMMARY 

We have analyzed the publicly available VL2 halo cata¬ 
logue in order to characterize the infall properties, orbital dy¬ 
namics, and structural evolution of over 6000 subhalos within 
a galactic host. Our main focus is on the role of subhalo con¬ 
centration in each of these categories, and how this relates to 
z = 0 observables, such as circular velocity and size. 

We define subhalo infall as the time when a halo reaches 
maximum mass. In other words, a halo becomes a subhalo 
when its growth is halted, mainly due to tidal truncation from 
the host. After infall, subhalos experience mass loss from tidal 
forces and exhibit internal readjustment as they gradually sink 
toward the host center. We focus on low-mass subhalos for 
which dynamical friction plays only a minor role, leading to 
qualitatively different behaviour than is often described for 
high-mass subhalos. In the following paragraphs we report 
the main results of our paper. 


Subhalo mass function: We compare the unevolved and 
evolved subhalo mass functions with the results published 
from the Aquarius simulations. The unevolved mass function 
uses the mass of each subhalo when it fell in while the evolved 
mass function uses the mass at z = 0 , showing the cumulative 
effect of tidal mass loss after infall. While the unevolved mass 
functions agree well, the evolved VL2 mass function is sys¬ 
tematically lower, by 30%, corresponding to a downward shift 
in mass (Figure!^. The lower normalization of the VL2 sim¬ 
ulation (cTg = 0.74 vs. erg = 0.9) could be the o rigin , although 
the physical explanation remains unclear (see §3.2|l. 


Properties at infall: Several important relationships 
among subhalo properties at the time of infall emerge; ( 1 ) 
The typical infalling halo mass does not evolve significantly 
with time, with a value of ~ 10^ Mq (Figure^ and also Fig¬ 
ure [Tfil l. ( 2 ) Rare, massive halos fall in muchlater than less 
massive ones, as expected in hierarchical structure formation 
(Figure]^. (3) Halos that fall in earlier have lower concen¬ 
trations, consistent with the well-known concentration-mass- 
redshift relationship for dark matter halos at fixed mass (e.g., 
[IGypin et al.||2011 see Figure 1^. (4) Halos that fall in ear- 
lier or, equivalently, that have lower concentrations at infall, 
tend to experience tidal growth truncation at a larger radius 
(Figures]^ and 1^. 


Energy and angular momentum: The orbital energy and 
angular momentum of subhalos at infall are significantly dif¬ 
ferent than reported in previous studies that focused on more 
massive subhalos. The low-mass subhalos found in VL2 are 
skewed toward lower specific binding energy (Figure [7ji and 
slightly more circular o rbits (Figure [9|l. This is consistent 
with the recent results of |Jiang et al.| ( |^15 1. The explanation 
for these trends may be related to the environment in which 
objects form - massive subhalos fall into the host preferen¬ 
tially along filaments, leading to tightly bound, radial orbits - 
and/or gravitational interactions prior to infall that preferen¬ 
tially inject energy and tangential motion into low-mass sub¬ 
halo orbits. 


Mass loss: We find that subhalos undergo most of their 
mass loss on the first pericenter passage, with a median mass 
loss fraction of ~ 0.2 (Figure [T0)i. The fraction of mass lost 
in the first pericenter passage is most correlated with halo 
concentration at infall - less concentrated halos tend to un¬ 
dergo more mass loss, nearly independent of mass (Figure 
[TT] l. There is a significant but less pronounced correlation 
oTmass loss fraction with subhalo mass, but this trend is only 
apparent, being explained by the fact that more massive ha¬ 
los have lower concentrations on average and are thus more 
susceptible to tidal effects. Scatter in mass loss at fixed con¬ 
centration can be mainly attributed to eccentricity with radial 
orbits tending to lose more mass than circular orbits. 

Orbital period: The period of the first orbit after infall is 
roughly proportional to the dynamical time of the host halo, 
foi-b oc Tdyn oc [A(z)pcrit(z)]~'^^- There is significant scatter, 
however, in the apocenter, rapo/Rvii, which results in a compa¬ 
rable scatter in forb/Tdyn (Figure [I^. The scatter originates in 
the concentration of the infalling subhalos; low-concentration 
subhalos begin to be disrupted earlier and thus experience 
much longer initial orbits than subhalos with higher concen¬ 
trations. 


Spherical Symmetry: Motion in a spherical potential, in 
which the direction of angular momentum does not change. 




















18 


is not a good approximation to subhalo orbital dynamics. In 
particular, the direction of the angular momentum vector is 
not fixed. After a few orbits, the direction of the angular mo¬ 
mentum is essentially randomized. This seems to be a generic 
feature of subhalo evolution in highly inhomogeneous, triax- 
ial host halos (Figure [T3]l. 

Evolution in the Vmax-^max plane: As subhalos are tidally 
disrupted by the host halo, their maximum circular velocities 
and radii steadily decrease, tracing out tracks in the Vmax-^max 
plane. While the joint median evolution, as well as their indi¬ 
vidual dependence on tidal mass, are in qualitative agreement 
with previous studies, we find a substantial amount of scat¬ 
ter. Furthermore, this scatter can be mostly attributed to vari¬ 
ations in the concentration at infall. The difference is most 
pronounced in the evolution of rmax^ subhalos that are more 
concentrated at infall experience a weaker evolution in r^ax 
as they lose mass (Figure [l4li. Concentration at infall deter¬ 
mines evolution in the Vmax-^max plane (Figure [T5]l. 

Disrupted Subhalo Population: We find that subhalos 
with early infall and/or low concentration are preferentially 
disrupted within the host (Figure [T6| ). We find no depen¬ 
dence of survivability on mass witEin the range probed by 
VL2. These results are consistent with the notion that tidal 
mass loss is correlated with concentration instead of mass for 
low-mass subhalos. The typical infalling mass of m ~ 10^ Mq 
( fiQ ~ 10 “^) is true for both surviving and disrupted subhalos. 
Circularity does not influence subhalo survivability. There is 
a slight bias in tightly bound orbits with p < 1 being prefer¬ 
entially disrupted though this is a relatively small effect. The 
inf all dist ribut ions an d fitting functions for p and e presented 
in |3.3.3|and ^3.3.4|based on surviving VL2 subhalos should 
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do a good job at representing the infall distribution for all sub¬ 
halos (surviving plus disrupted) that ever fell onto the host 
(Figures and [TSl l. This is the case at both Zinfaii with an 
NFW host potential as well as Zcross with a point mass host 
potential. 

In this work we have extended previous detailed analyses 
of subhalo dynamics and evolution to the much lower mass 
ratios probed by the Via Lactea 11 data. We have found quali¬ 
tatively different behaviour in this low-mass regime, with dy¬ 
namical friction and orbital dynamics playing a lesser role, 
and the interior structure of the subhalos, expressed in terms 
of concentration, playing a much more important role. 

The ‘concentration bias’ we And here raises the prospects 
of significantly improving our ability to connect ultra-faint 
dwarf galaxies to the primordial fluctuations from which they 
collapsed. More detailed study, in particular with finer time 
resolution and a larger sample of simulated Galactic host ha¬ 
los, will be necessary before we can reliably use concentration 
bias in near held cosmology. 
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